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Abstract 

In this paper we propose a method to study quantitatively the glass transition in 
a system of interacting particles. In spite of the absence of any quenched disorder, 
we introduce a replicated version of the hypernetted chain equations. The solution 
of these equations, for hard or soft spheres, signals a transition to the glass phase. 
However the predicted value of the energy and specific heat in the glass phase are 
wrong, calling for an improvement of this method. 
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1 Introduction 



The long lasting studies of the glass transition, in spite of recent progress |I| , have not yet 
produced a completely convincing theory deriving the low temperature behavior from a 
microscopic Hamiltonian of interacting particles. 

It has been proposed that near this transition real glasses could behave in the same 
way as some disordered systems which are described by a Hamiltonian with quenched 
randomness ^, |, . This proposal is supported the existence of a wide class of mean 
field systems, with fixed non random Hamiltonians which show a typical glassy behaviour 

i> B B HO' O' 0' mi- Although they depend on one or a few coupling constants. 



these systems behave as if they were typical samples of a class of systems with quenched 
randomness (apart from the possible existence of a crystal phase 

The relevant disordered systems found so far have the following behaviour at the mean 
field level. There are two phase transitions: a usual static transition at low temperature 
(Tr), where the replica symmetry is broken and the specific heat is discontinuous, and a 
dynamic transition at a higher temperature (Td), where some time persistent correlations 
set in. In the region where T < Td the infinite system, quenched from some high energy 
configuration, gets trapped into metastable states, such that its energy density is larger 
than the equilibrium value. In this dynamics the metastability appears at Tr and the 



thermodynamic transition, which is present at a lower temperature, is not accessible 0, 
p!5| , |16|, |18[. Beyond the mean field approximation one finds that the equilibration time is 
not divergent at To, but it becomes exponentially large in the region Tq < T < T^. It has 
been argued that in some spin model the equilibration time diverges (as exp{A\T — Tc\~'^\ 
7 = (i — 1 in (i dimensions 0, ^ . This problem will not be studied here, we shall stay at a 
mean field level. 

In this paper we want to present a progress report on an analytical study of the glass 
transition for a real three dimensional model of interacting particles. 



2 Breaking the replica symmetry 

We consider a system of interacting particle in a volume V. We study the infinite volume 
limit in which ^ oo at fixed p = N/V. The Hamiltonian is given by: 

H{x) = J2U{x^-Xk) (1) 

i<k 

We will consider two forms of the potential among particles: The hard sphere case where 
U{x) = oo for r = \x\ < 1 and U{x) =0 for r > 1, and a soft sphere case U{x) = r~^^. 

In both cases, if the system is cooled from high to low temperature fast enough, crys- 
talization is inhibited and the system undergoes a glass transition. For hard spheres this 
transition occurs at a density of about 1.15 independently from the temperature, while 
the crystallisation would occur for an equilibrium system at density around .95. In the 
case of soft spheres thermodynamic quantities are functions of the dimensionless density 
7 = pT^~^/^^ and the glass transition is located around 7 ~ 1.6, while the freezing value of 



7 is around 7 ~ 1.15 [IS 



It is possible to consider n replicas of the same system and to study the total partition 
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function 

Z^JUdxt exp -P H{x^) + e ^ ^ - x\)\ . (2) 

i,a \ a. a,b i^k ) 

The detailed form of the (attractive) potential ^a,h is not relevant here because we will 
eventually set e = 0. 

Replica symmetry breaking corresponds to the case in which (in the infinite system) 
the correlation among replicas does not vanish in the limit e — > 0. The physical meaning 
of this breaking is clear: at high density, if several metastable states contribute to the 
Boltzmann measure, different replicas may be frozen in the same state. Indeed we can 
propose the following picture of the glass phase. 

There are many equilibrium states, which we label by an index a (modulo rotations and 
translations). These states are identified by the density < p{x) >a= Pa{x) which depends 
on X in each state. In this paper we denote by p{x) the density operator X)j ^{^ ~ Xj) 
and similarly by Pa{x) the density operator in replica a. The brackets stand for thermal 
expectation values. Each state a appears in the Boltzmann measure with a weight Wa- We 
assume that there are many states, and the quantities which we can access are averaged 
over the states. Since the one point average J2a '^aPa{^) = p is uniform, we need to study 
the two point correlations like: 

9{x,y) = \^Woc^< 5{xi- x)5{xj - y) >oc 

gi{x,y) = -Y^^IY, < 6{xi- x) >a< 6{xj -y) >a , (3) 

As for the densities in two different states, a ^ (3, we assume that they are uncorrelated. 

In replica space one could thus expect the correlations to be, for a ^ b: < Pa{x)pb{y) >— 
'W2Pa{x)Pa{y) — 'W2P^gi{x,y), wherc we use the notation 

Wk =< Xl'^a / P{w) w^. (4) 

Notice that we are describing the quantities w by their probability distribution P{w). This 
description may look strange for such a deterministic problem (the Hamiltonian is not 
random). It is clear that in a finite volume, at fixed N, the w's are given numbers. However 
they may strongly depend on N so that if we average over a small window in N around N 
and send — > oo we induce a smooth probability distribution for these variables. In this 
formulation one finds the unpleasant feature that the correlation function of more replicas 
do not factorize when the distance goes infinity. For instance, the four point function 

< pa{x)pb{y)pc{z)pd{t) >= ^W^Pa(x)pa(y)Pa(^)Pa(^) (5) 

a 

goes to W4gi{x,y)gi{z,t) when x — y = 0(1), z — t — 0{1) and x — z ^ oo, which is 
different from the factorized form {w2ygi{x,y)gi{z,t). 

The non factorization of the correlation function is deeply linked to the probability 
distribution of the variables w. This crucial ingredient of the description is usually dealt 
with by using the replica symmetry breaking formalism. This means that the correlation 
functions (e.g. < Pa{x)pb{y) >) are no more symmetric under the permutation of the 
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replicas, but they have the advantage of being clustering [|. This formalism provides and 
automatic bookkeeping of all complications which would arise from the existence of many 
states. If one assumes a given form for replica symmetry breaking, it corresponds to a 
given form for the probability distribution of the w. In the simplest case, called one step 
breaking, one divides the n replicas in n/m groups of m replicas and one assumes the 
following structure of correlation functions: 

< pa{x)pa{y) > = p^g{x - y) , 

< pa{x)pb{y) > = p^gi{x — for a 7^ 6 in the same group , 

< pa{x)pb{y) > = p^ for a and b in different groups. . (6) 

In the physically relevant case where m < 1, this form of the correlations corresponds 
to assume that Wa = Gxp{—i3fa)/[J2a'^w{~Pfa')] where the fa are negative quantities 
extracted with a probability distribution j3mexp{l3mf). The details of this construction 
are described extensively in the literature [^. 

May be the reader has some familiarity with the replica formalism and he is waiting here 
for the canonical sentence we send n to zero at the end of the computation. However this is 
not a random system and there is nothing to average (apart from the value of N), there is 
no need to send n to zero. Indeed it turns out that in this scheme, the free energy density, 
— \ia{Z)/{i3Nn), is independent of n as soon as < pa{x)pb{y) > — < Pa{x) >< phiy) >= 
for a and b in different groups. In some sense n is infinite, because we assume that m 
divides n for any m. 

In this scheme the free energy depends on m, which should be considered as a variational 
parameter and the free energy must be maximized (not minimized!) as usual in the replica 
approach. For m = 1 we recover the replica symmetric result. The static glass transition 
temperature Tr is characterized by the existence of a maximum of the free energy at 
m(T) < 1, where ^(Tr) = 1. On the other hand, the appearance of metastable states 
which will trap the dynamical evolution of the system is signalled by the existence of a non 
trivial gi{x) in the limit m ^ 1 . 



The programme is thus clear. The real difficulty consists now in implementing it, i.e. 
in computing the properties of our replicated system of interacting particles, allowing for 
replica symmetry breaking. A first step in this direction will be proposed in the next 
section. 



3 The replicated hypernetted chain approximation 

The hypernetted chain (HNC) approximation consists in considering only a given class of 



diagrams in the virial expansion It gives a reasonable account of the liquid phase, 

and it has also been used for studying the first order transition to the crystal phase P^ . 
We will consider here this approximation because it has the advantage of having a simple 
variational formulation. In the liquid phase, where the density is constant, the usual HNC 
equation (for the non replicated system) can be written as 

g{x) = exp {-f3U{x) + W{x)) , (7) 



^We thank Francesco Guerra for discussions on this point. 
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where: 



p'g{x) = p\l + h{x)) =< p{x)p{Q) > -p5{x), 
Wix) = [J^e-^^^^^^ (8) 

and we denote by h(p) the Fourier transform of h{x). 

This equation can be derived by minimizing with respect to h{x) the following free 
energy per unit volume, in the space of functions of 

2/3F = J d'x p'g{x)[Hg{x)) - 1 + pU{x)] + j ^L(ph(g)), (9) 

where L{x) = — ln(l + x) + x — x'^/2. 

In this note we propose a bold generalization of the HNC equations for n replicas. The 
replicated free energy is now 

2/3F = p'J d''xY,gab{x)[ln{g,b{x)) - 1 - f3{Uix)Sa,b + e<Pab)] + TrL(ph), (10) 

a,b 

where h is now an operator both in x space and in replica space. 

If in the limit e one finds that at large enough densities gab is non zero off the 
diagonal, replica symmetry is broken. In the case where gab is of the form shown in eq. (^), 
this equation can be used to compute the properties of the correlation function in the glassy 
phase. This approach amounts to a study of the density modulations in the glass phase 
at the level of the two point function. In the glass phase pa{x) becomes space dependent. 
However, as argued in the introduction, the necessity of averaging over the states a forces 
us to study this x dependence at the level of the two point correlations: We keep here to 
the definition of g through J^a'Wa < p{^)piy) >a= S{x — y) Y^a'>^aPa{x) + p'^g{x,y), while 
the second term should be written more properly as pa{x)ga{x, y)pa{y)- So our correlation 
g reflects the structure of pa{x) as a sum of peaks of unit weights, smoothed by the average 
over states. 

Within the one step breaking scheme (|^), the free energy is: 

hm ^ = / d''x{~gix)[H~gix)) - 1 + /3f/(x)] - (1 - m)g,ix)[\nig,ix)) - 1]} 
np J 



{27tY \ 2 

1 ~ 1 — TTl — \ 

+^ ln[l + phiq) - (1 - m)phiq)] T ^^l^ + P^(^) " (^1) 

The static transition is identified as the temperature (or density) at which there exists 
a non trivial solution to the replicated HNC equations dFjdgix) = 0, dF/dgi{x) = and 
dF/dm = 0, for m G [0, 1] (in fact we must minimize the free energy with respect to g, but 
maximize with respect to gi and m). The equations to be solved for the statics are thus: 

g{x) = exp (^-i3U{x) + W{x)'j , gi{x) = exp {Wi{x)) 

d'x g,{x)[Hg,{x)) - 1] = f + ^ (12) 

^ ^^^ ^^ ^'^ ^ J i2nY\ 2 mpl+ph(g)-(l-m)phi(g) ^ ^ 

1 l + ph(g)-(l-m)phi(g)\ 
;272 In . . s I (13) 



1 + ph(g) -phi(g) 
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where the Fourier transforms W(g) and Wi(g) of W{x) and Wi{x) satisfy: 

(h(g)-(l-m)hi(g))2 



W(g) - (l-m)Wi(g) = p 



l + p(h(g)-(l-m)hi(g)) 



^ . 1 / p(h(g)-(l-m)hi(g))^ _ p(h(g) - hi(g))^ \ 

m\l + p{h{q)~{l-m)h^{q)) 1 + p(h(g) - hi(g)) 7 

The dynamical transition may be characterized as the highest temperature at which 
there is a non trivial solution of the stationarity equations dF/dg{x) = and dF/dgi{x) = 
at m = 1~. The corresponding equations are obtained by substituting m ^ 1 in the 
first two equations of (]T3|). The equation for g is identical to the usual HNC equation (|), 
while gi is a solution of gi{x) = exp(iyi(a;)), with Wi given by the second equation of ( p!4D 
at m = 1. 

Within this static approach there is no obvious definition of the dynamic energy. Previ- 
ous work ^] suggests the following computation: one considers a system at equilibrium 
at a temperature Tr > Tq and an other system at temperature T. One defines the corre- 
lation functions 

P^9r{x -y) = < p{x)p{y) >r, 
p'^gix -y) = < p{x)p{y) >, 
p'^9o{x -y) = < p{x) >a (< p{y) >ij)„, 

P^gii.x-y) = < p{x) >a< piy) >a, (15) 

where we have assumed that the states of the system at a given temperature are in corre- 
spondence with the states at lower temperature. One imposes the constraint on the system 
at temperature T that the function gQ and gi are not zero. If T = T^j this constraint can be 
easily satisfied. If T < Tr this constraint forces the system at T to stay out of equilibrium. 



It has been conjectured in |16| that the correlation functions in the metastable state (from 



which we can easily extract the energy) are given by gi{x) computed at Tr = Tm{T). 



4 Solution of the equations 

Let us first discuss some technical points which are common to the soft and the hard sphere 
cases. In both cases the first task is to solve the replica symmetric HNC equation. For 
spherically symmetric functions in dimension three we use the Fourier transform for the 
radial dependance, in the following form: 

POO 

gh(g) = 27r / dr sm{qr)rh{r). (16) 
Jo 

We discretize this formula introducing in r space a cutoff R and a mesh size a. In this 
way we have a simple formula for the inverse Fourier transform and we can also use the fast 
Fourier transform algorithm. In most of the computations we have taken a = 1/32.5 and 
L = 128 * a ~ 4. We have tried smaller values of a and larger values of L without serious 
effects. The solution of the equations can be found either by using a library minimization 
program 0, or a program which solves non linear equations. We have found first the solution 

^When one has to minimize the free energy with respect to one variable and maximize it with respect 
to another one, we first minimize with respect to the first variable, later maximize with respect to the 
second, and iterate the procedure until convergence. 
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at low enough density and then followed it by continuity while gradually increasing the 
density. 

4.1 The soft sphere case 

The HNC equation gives a description of the liquid phase which is not perfect, but pre- 
cise enough for our purpose: The energy (or equivalently the pressure), does not depart 
more than 15% from the correct value (see fig.l), and the correlation function is also well 
reproduced (see fig. 2). 

The numerical solution of the replicated HNC equations finds a dynamical transition 
at 7 ~ 2.05, and a static replica symmetry breaking solution at 7 ~ 2.15. In numerical 
simulations the glass transition is found at a smaller value of 7, namely 7 = 1.6. In the glass 
phase, the correlation function gi{r) is essentially a smoothed form of the function g{r) 
plus an extra contribution at short distance which has integral near to 1 (see Fig.3). This 
form seems very reasonable: Considering the definition (^) of gi, we see that it basically 
characterizes the average over a of the product pa{x)pa{y), which is precisely expected to 
have this kind of peak structure. 

In spite of this nice form for gi, this solution has some problems. A first one is found on 
the value of the energy. The static energy as function of 7 (or equivalently of T at density 
1). is plotted in fig.l. Although there is a discontinuity in the specific heat, it is extremely 
small and the final effects on the internal energy are more or less invisible. The specific 
heat remains extremely large. Moreover the value of m has a very unusual dependence 
on the temperature (see fig.4). In all the known models with one step replica symmetry 
breaking, the breakpoint m varies linearly with T at low temperatures. Here we have a 
very different behaviour. We have also computed the dynamical internal energy and found 
out that it differs from the equilibrium one by an extremely small account. 

We conclude that if we consider the qualitative behavior of the correlation functions, 
we find a reasonable form, on the other hand the energy in the glassy phase turns out to 
be quite wrong. 

4.2 The hard sphere case 

The situation is quite similar in the case of hard sphere. The HNC approximation works 
reasonably in the liquid phase. The pressure, which can be extracted either from the free 
energy or from its relation with g{l~^), does not depart from the correct one by more than 
15% in the liquid phase. The main defect is the absence of the peak around r = and 
a too large value of the peak at r = 2. The pressure seems to diverge proportionally to 
(Pc — p)^) where pc is around 1.6, while the maximum possible density, corresponding to 
the fee lattice, is \/2. 

Here replica symmetry is broken around p ^ 1.19 and the dynamical transition is 
located at p ~ 1.17. These values are very close to the result of the numerical simulations 
which find a freezing transition around p ~ 1.15. Unfortunately also in this case the 
computed value for the pressure differs by a very small amount from the replica symmetric 
one and it is therefore unacceptable. Similar conclusions have also been reached for the 
Lennard Jones potential. 
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5 Conclusions 



The simple implementation of the replica approach to glasses which we have proposed here 
provides some interesting results, like the existence of a glass transition at a reasonable 
value of the density. However it is not satisfactory, in the sense that the effects of the 
transition on the thermodynamic of the systems are much too small. 

At the moment we do not have a clear understanding of the reasons of this failure to 
grasp the thermodynamic properties of the glassy phase. Two possible explanations came 
to our mind. The first one would be that an approximation like that of HNC may miss 
some of the main physical characteristics of the problem in the glassy phase. A first look at 
the corrections to the free energy indicates that they are quite large, lending some support 
to this hypothesis. More work is needed to decide on what class of diagrams should be 
added to cure this problem. On the other hand we must also admit that we have not found 
so far a full derivation of the replicated HNC equations. A more proper approach would 
be to work with a free energy expressed as a functional of the density p{x), seeking all the 
(non translational invariant) solutions corresponding to glass phases. Then the correlations 
and thermodynamics could be studied by giving to each solution a weight proportional to 
exp(— Fa being the free energy of the solution labeled by a. 

We think that finding an analytic approximation scheme which produces reasonable 
results in the glass phase is within reach. The method we propose seems to be promising 
in this respect and shows how the replica method could be used to study the glass phase. 
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Figure Captions 

Fig.l The static energy of a system of soft spheres as a function of the dimensionless inverse 
density 7 from numerical simulation (squares), from the replica symmetric HNC equation 
(full line) and from one step replica symmetry broken HNC equation (open circles). 
Fig.2 The correlation function g{r) of a system of soft spheres as a function of the distance 
at the dimensionless density 7 = 1.6 corresponding to the numerical glassy transition: 
numerical simulations (points) and replica symmetric HNC equation (full line). 
Fig. 3 The correlations g{r) (full line) and gi{r) (broken line) obtained from the replicated 
hnc equations, as functions of the distance, for soft spheres at the density where replica 
symmetry is broken (7 = 2.15). 

Fig.4 The value of the breakpoint m in the replica symmetry breaking solution for soft 
spheres as a function of the dimensionless density 7. 
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